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We consider the motion of classical spinning test particles in Schwarzschild and Kerr metrics and 
investigate innermost stable circular orbits (ISCO). The main goal of this work is to find analytically 
the small-spin corrections for the parameters of ISCO (radius, total angular momentum, energy, 
orbital angular frequency) of spinning test particles in the case of vectors of black hole spin, particle 
spin and orbital angular momentum being collinear to each other. We analytically derive the small- 
spin linear corrections for arbitrary Kerr parameter a. The cases of Schwarzschild, slowly rotating 
and extreme Kerr black hole are considered in details. For a slowly rotating black hole the ISCO 
parameters are obtained up to quadratic in a and particle’s spin s terms. From the formulae 
obtained it is seen that the spin-orbital coupling has attractive character when spin and angular 
momentum are parallel and repulsive when they are antiparallel. For the case of the extreme Kerr 
black hole with co-rotating particle we succeed to find the exact analytical solution for the limiting 
ISCO parameters for arbitrary spin. It has been shown that the limiting values of ISCO radius and 
frequency do not depend on the particle’s spin while values of energy and total angular momentum 
depend on it. We have also considered circular orbits of arbitrary radius and have found small-spin 
linear corrections for the total angular momentum and energy at given radius. System of equations 
for numerical calculation of ISCO parameters for arbitrary a and s is also explicitly written. 

PACS numbers: 04.20.-q - 04.25.-g 


I. INTRODUCTION 

In General Relativity the presence of a rotation (intrin¬ 
sic angular momentum, or spin) of the central body in¬ 
fluences motion of a particle orbiting it. Due to this rea¬ 
son the orbits of test particles differ in the Schwarzschild 
and Kerr backgrounds. When, in turn, a test particle 
has spin as well, it will also influence the particle’s orbit. 
In particular, the motion of a spinning particle will dif¬ 
fer from the non-spinning one even in the Schwarzschild 
background. 

Let us consider a massive non-spinning test particle 
orbiting a central black hole (BH). There exists a mini¬ 
mal radius at which stable circular motion is still possi¬ 
ble, it defines the so-called innermost stable circular or¬ 
bit (ISGO) in given background. For the Schwarzschild 
background the radius of ISCO equals to 6iwQ dll. In 
the Kerr space-time circular motion is possible only in 
the equatorial plane of BH and the radius of ISCO de¬ 
pends on the direction of motion of the particle in com¬ 
parison with the direction of BH rotation, whether they 
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which will be introduced further have the following dimensional¬ 
ities: [L] = [M], [J] = [M], [E] = 1, [a] = [M], [s] = [M], 


co-rotate or counter-rotate. Co-rotation and counter¬ 
rotation cases correspond to parallel and antiparallel ori¬ 
entation of vectors of the orbital angular momentum of 
the particle and the BH angular momentum. For the case 
of the extreme Kerr background the difference between 
these two variants is quite considerable: we have 9M for 
the antiparallel and M for the parallel orientation 

Values of the ISCO parameters (radius, total angular 
momentum, energy, orbital angular frequency) are deter¬ 
mined by the Kerr parameter a. The main subject of the 
present paper is to investigate how the ISCO parameters 
(at a given value of a) are changed if a test particle has 
spin (Fig. [T|). We consider motion of a spinning test 
particles in the Schwarzschild and Kerr metrics in the 
equatorial plane and analytically derive the corrections 
for the parameters of ISCO taking spin to be a small pa¬ 
rameter. We restrict our consideration by the case when 
vectors of black hole spin, particle spin and orbital angu¬ 
lar momentum are collinear. 

For non-spinning particles moving in the Schwarzschild 
metric the ISCO parameters were found by Kaplan [I|. 
The solution of Einstein field equations around a rotat¬ 
ing black hole was found by Kerr Q. In the following 
works of Carter @ and de Felice Q the geodesic mo¬ 
tion of general type was studied for the Kerr and Kerr- 
Newman (charged BH) geometry. See also Wilkins Q 
for bound orbits in Kerr metric and Dymnikova @ for 
review. The parameters of ISCO in Kerr space-time for 
a non-spinning particle were obtained in the works by 
Rufhni & Wheeler Q and Bardeen [3|. 
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The problem of the motion of a classical spinning 
test body in General Relativity was considered in pa¬ 
pers of Mathisson 0, Papapetrou m and Dixon [l^ 
Il4| . using different techniques. The equations of motion 
of a spinning test particle in a given gravitational field 
were derived in different forms; they are now referred to 
as Mathisson-Papapetrou-Dixon equations. From these 
equations it follows that the motion of the centre of mass 
and the particle rotation are connected to each other, 
and when the particle has spin the orbits will differ from 
geodesics of a spinless massive particle. In this paper 
we use the equations of motion in the form derived in 
Dixon’s papers. 

Influence of spin on orbits in the Schwarzschild metric 
was investigated in the paper of Corinaldesi and Papa¬ 
petrou [l^ , and in the paper of Micoulaut . A motion 
of a spinning test particle in Kerr metric was considered 
by Rasband El. Using integrals of motion arising from 
the symmetries in the Kerr space-time Rasband [l^ de¬ 
rived the equations of motion of the radial coordinate 
of a spinning particle for the motion in the equatorial 
plane of a BH. Also in that work the spin-induced correc¬ 
tions for the radii of last stable orbits for Schwarzschild 
(a = 0) and extreme Kerr (a = M) metric were illus¬ 
trated through numerical calculations. In the paper by 
Tod, de Felice & Calvani [l^ the influence of test body’s 
spin on the radii of ISCO for different values of Kerr pa¬ 
rameter a including naked singularities was numerically 
calculated, see also papers of Abramowicz & Calvani 
and Calvani More general case of Kerr-Newman 

metric was analysed by Hojman & Hojman [2l|. Subse¬ 
quently the number of works on this subject were pub¬ 
lished (23 - 14^ . The detailed derivation of the equations 
of motion using the integrals of motion was presented in 
the work of Suzuki & Maeda for the Schwarzschild 
case and in the work of Saiio et al. [2^ for the Kerr back¬ 
ground. Tanaka et al. suggested that the radius of 
the last stable orbit is independent of the particle’s spin 
in the extreme Kerr background (for co-rotating orbits). 

Many authors used an effective radial potential for in- 
vesti gati on of spinning particle orbits [ 13 , [Il|, [2ll - l^ . 
[^. jdljl. [13 . Method of calculation of ISCO param¬ 
eters of spinning particle moving in Kerr metric is pre¬ 
sented in details in the paper (39| . The equations for 
the circular orbits in the equatorial plane resulting from 
the Mathisson-Papapetrou-Dixon equations appear to be 
irresolvable analytically, and ISCO parameters for arbi¬ 
trary value of a are supposed to be found numerically. 
Linear corrections in spin for the ISCO parameters in 
Schwarzschild metric were found by Favata [s^. In this 
work we analytically obtain the small spin corrections 
for the ISCO parameters for the Kerr metric at arbitrary 
value of a. 

For the extreme a = M and almost extreme a = 
(1 — S)M Kerr BH we succeed to find the exact ana¬ 
lytical solution for the ISCO parameters for arbitrary 
spin, with only restrictions connected with applicability 
of Mathisson-Papapetrou-Dixon equations. It has been 
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FIG. 1. The influence of BH and test-body’s spin on the radii 
of ISCO. 

In case of Schwarzschild BH the ISCO radius for spinless par¬ 
ticles equals to 6M, see thick solid line. For a particle with 
spin, ISCO radius splits in two cases, with different mutual 
orientation of spin and orbital angular momentum (parallel 
or antiparallel), see dash lines. 

In the case of Kerr BH and spinless particle, ISCO radius 
splits in two cases, with co- and counter-rotation to BH (an¬ 
gular momentum is parallel or antiparallel to BH spin), see 
thin solid lines. For a particle with spin each line additionally 
splits in two dashed lines corresponding to different spin ori¬ 
entation. The only exception is the case of extreme Kerr BH 
with a co-rotating particle: in this case the ISCO radius is 
equal to M and does not depend on magnitude and direction 
of particle spin. 


shown that the limiting values of ISCO radius and fre¬ 
quency for a = M do not depend on the particle’s spin 
while values of energy and total angular momentum do 
depend on it. 

The present paper is organised as follows. In section 
m we describe a motion of a test body without spin and 
introduce the notion of the effective potential (EP). In 
section lllll we describe how to find ISCO parameters with 
the help of EPs, in the Schwarzschild and Kerr metric. In 
section m we present the basic equations for motion of 
a spinning test body. In section |V] we introduce EPs for 
spinning particles and write explicitly equations for ISCO 
parameters for arbitrary Kerr parameter and arbitrary 
spin. In section EH we derive formulae for calculation of 
the small-spin linear corrections to ISCO parameters for 
arbitrary value of a. In sections fVIIl I VIlIl IIXI we consider 
in details the cases of Schwarzschild metric, Kerr metric 
with a M and extreme Kerr metric. In section we 
present an exact solution for extreme Kerr BH in the case 
of co-rotation. In section m we consider circular orbits 
of arbitrary radius. Section IXIII is conclusions. 
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II. MOTION OF A SPINLESS TEST BODY. 
EFFECTIVE POTENTIAL 


In the problems of motion of a test particle in a central 
field it is convenient to express the equations of motion 
for the radial coordinate in terms of a so-called ” effective 
potential” energy. 

In the Newtonian dynamics these equations for a par¬ 
ticle moving in the gravitational field of a central body 
with mass M can be written as (see, for example, 

K9 


where e is the total non-relativistic energy per unit mass, 
C/tv is the so-called effective potential (EP) 


Un = 



( 2 ) 


Here L is the angular momentum per unit mass, G is the 
gravitational constant. 

In the framework of General Relativity, the 
Schwarzschild space-time is given by the expression 
for the interval (in spherical Schwarzschild coordinates 
as: 


ds^ 




+ {de^ + sin^ e dcp^) . 


( 3 ) 


energy. Newtonian non-relativistic energy and effective 
potential are obtained as e = (£1^ — l)/2 and Un = 
(Uschw - l)/2, see [ii|. 

Let us consider the effective potential in the 
Schwarzschild space-time © as a function of radial coor¬ 
dinate. Then, the angular momentum L plays a role of a 
parameter defining the shape of Uschwir) curve. In the 
case L > 2y/3M the effective potential has two extrema: 
maximum and minimum, at the radii of which unstable 
and stable circular motion are possible correspondingly. 
In the case when L equals to the boundary value 2y/3M, 
two extrema of EP merge into one inflection point. This 
boundary value of L defines parameters of the last stable 
orbit which is also called the innermost stable circular 
orbit (ISCO), i. e. the boundary orbit on which the fi¬ 
nite motion is still possible. Further, when the angular 
momentum L of a test body is less than 2^/3M, the EP 
does not have an extremum. For these values of angular 
momentum neither type of finite motion is possible and 
a test body will inevitably fall in the black hole whatever 
values of E it may have. 

In Boyer-Lindquist coordinates the Kerr metric is given 
by the expression [3, l45j| 


ds^ = -{l- 


2Mr\ , 9 AMarsin^d , , E , , 


dt — 


-dt dip + -^dr 
A 


^ ,/,2 / 2 2 2 Mra?'s\r? 9 \ .2^,2 

+ E d 9 ^ + { U + -\ -- I sm^ 9 dip^, 


( 8 ) 


In this space-time the equation of motion for the radial 
coordinate is 



where a is the specific angular momentum of a black hole, 
S = -I- cos^ 0, A = — 2Mr + ■ 

The equation of motion for the radial coordinate of a 
test particle in the equatorial plane (9 = tt/2) of Kerr 
metric is written as [45l | 


with the connection between dt and dr as 

dt E 

dr 1 — 2M/r 


( 5 ) 


Here E is the total energy at infinity per unit particle 
rest mass, L is the angular momentum per unit particle 
rest mass, t is the proper time. This equation can also 
be presented in terms of an effective potential 


?:i =E'‘-u^ 


dr 


Schw * 


( 6 ) 


Here we identify with the effective potential the following 
expression: 


Uschw 



( 7 ) 


Comparing equations 0 and (j4]) we should keep in mind 
that non-relativistic total energy e is not an analogue of 
E, since unlike the latter, it does not contain rest mass 



= 


2M a^(£;^-I)-L^ 


2M{L-aEf' 


(9) 


with the connection between dt and dr 


dt 

dr 


2Ma 

rA 




2Ma^ 


( 10 ) 


Here L = is a z-component of the angular momentum 
per unit rest mass, equal to the total angular momentum 
per unit mass, when z-axis is parallel to the axis of the 
rotation of BH. So L > 0 and L < 0 corresponds to the 
co-rotation and the counter-rotation cases respectively, 
and a is always positive. 

We see at the right-hand side of the equation (0 that 
the expression in brackets is dependent on the energy E. 
For this reason the effective potential in the Kerr space- 
time cannot be defined in such a simple way as for the 
Schwarzschild space-time. 
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Let us write the equation for the radial motion as 


what corresponds to the equality 


r 


4 



{r^ + + 2Mra^)E‘^ - 


- AMrLaE + 2ML^r - + 2Mr^. 

( 11 ) 


It is convenient to write this equation as 


HI 


LI — USchw ■ 


(17) 


(ii) the acceleration of the radial coordinate should be 
absent: 



(18) 


J) =^[aE^-2PE + ^), 


( 12 ) 


Differentiating (|6]) with respect to r and dividing by f, 
we obtain 


where a, /3 and 7 are defined as: 

a = (r^ + - Aa^ > 0, 

^ = [o (r^ + a^) — Aa] L = 2MraL , 
7 = - A(r 2 + L") . 


Let us rewrite dED as 


= a{ E- 


P + - aj 


X \ E- 




a7 


and define EP equal to 


(13) 


> (19) 

so the condition (fTSl) corresponds to the relation 

^^[E^-Um=0. (20) 

For a given E we have 


(14) 


^^ = 0 , Usch^>0. ( 21 ) 

dr 

Solving the system of equations and (EU, we ob¬ 
tain the expressions for the radius of the circular orbit r 
and for the particle energy E on it , as a function 
of L in the form 


TT 7 -^ _ “ “'I 

UxerrytL/) — . (1^) 

Q 

The positive square root must be taken here, see 
( 4 ^ . This definition keeps analogy with the EP in the 
Schwarzschild space-time (and, therefore, with the New¬ 
tonian dynamics as well). 

In some other works, however, d, 113 the entire 
right-hand side of the equation (EJ is used instead of the 
effective potential Uxerr, especially for the investigation 
of circular orbits and the ISCO parameters. Despite 
that it contains energy, some authors d) IHl keep calling 
it an ’effective potential’, because it has some features of 
the EP important for the investigation of circular orbits 
(see next section). 


III. ISCO FOR NON-SPINNING PARTICLES IN 
THE SCHWARZSCHILD AND KERR 
SPACE-TIMES 


Let us describe how to find ISCO parameters and first 
consider the Schwarzschild metric. 

For a circular motion we need two conditions to be 
satisfied simultaneously: 

(i) the radial velocity should be equal to zero: 



(16) 


r 


2M 


1± 


I 12M2 


( 22 ) 


E = L 




(23) 


In eq. (1^ the upper sign corresponds to the stable cir¬ 
cular orbits in the minimum of EP and the lower one 
stands for the unstable circular orbits in its maximum. 
In eq. (l23l) radius r given by (1221) should be substituted. 
Prom the expression we see that for L < 2\/iM the 
radius becomes complex and we do not have a circular 
orbit at all. Mathematically this means that the EP with 
such values of L does not have extrema. 

In order to find the last stable orbit (ISCO) we need 
to find the last value of L at which the EP still has ex¬ 
tremum. The ISCO takes place when points of maximum 
and minimum of EP merge. Therefore we need the third 
condition: 

(iii) the inflection point of the EP satisfies relation: 


^2 tt 2 

“ Schw ^ n 


(24) 


So, in order to find the parameters of ISCO (r, E, L), we 
need to solve three equations (fT71) . (1^ . (IMl) simultane¬ 
ously. 
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From (I22p it is evident that the points of maximum 
and minimum of EP merge at the marginal angular mo¬ 
mentum L = 2'/3M. Thus, the ISCO parameters in the 
Schwarzschild space-time are [I|: 

’’ISCO = 6M, 

^isco = 2\/3M, 

Eisco = 

In the Kerr space-time the finding of ISCO parameters 
is more complicated. ISCO in the Kerr space-time were 
considered in i, i , where, in particular, parameters of 
ISCO for the extreme Kerr BH were found. This problem 
is described at length, for example, in the textbook by 
Hobson et al. [45j| . 

It is convenient to solve the problem by introducing a 
function 



Vir;L,E) = ^iaE^-2pE + j). (26) 

The function V{r;L^E) is the right-hand side of eq. 
(HH and has qualities which are important in our re¬ 
search. The expression for acceleration is given by the 
first derivative r = {l/2)dV/dr (compare with (fT^ l. and 
stability of the circular orbit is given by the sign of its 
second derivative d'^V/dr'^. It is convenient to use vari¬ 
able u = 1/r instead of r and variable x = L — aE instead 
of L. 

For the circular motion we need the conditions for the 
velocity (HU) and the acceleration (fT51) to be satisfied si¬ 
multaneously. Therefore for the Kerr metric we have 
from (|M|l the system of equations 


r v^ = o, 

< ^ _ 0 (27) 

I dr 

The usage of u instead of r does not change anything in 
the form of the system: since dV/dr = {dV/du){du/dr), 
the condition dV/dr = 0 is equivalent to dV/du = 0. 

The solution of the system (|27|l defines the set of pa¬ 
rameters X = L — aE and E for stable and unstable cir¬ 
cular orbits as functions of u [1^ . For the stable circular 
orbit we have 


x(u) = — 

E{u) = 


ly/u ± 


1 1/2 ’ 


u(l- 3Mu T 
1 - 2MuT ay/Mu^ 

Using L = a; -I- aE, we have: 

•\/M(1 -|- a?v? ± 2ay/Mu^) 


L(u) = T- 


u — 3Mzi =F 2a\/Mu^ 


1/2- 


(28) 


(29) 


The upper sign corresponds to the antiparallel orienta¬ 
tion of particle’s angular momentum L and BH spin a 
(counter-rotation), the lower - to the parallel one (coro¬ 
tation) . 

In order to find the parameters of ISCO, we need to 
add the condition 


d^V 

dr"^ 


= 0 . 


(30) 


Written in terms of u the third condition transforms to 



dr"^ dv? V dr / du dr'^ \ du"^ du 

(31) 

Since dV/du = 0 for a circular orbit, the third condition 
is d^V/du^ = 0. 

The three equations U = 0, dV/du = 0, d'^V/du^ = 0 
form a closed system on three parameters of ISCO E, x 
and u, which are then dependent only on M and the Kerr 
parameter a. 

Solving this system , we obtain the equation for the 
inverse ISCO radii in the Kerr metrics: 

I — 3a^M^ — 6uM =F 8ay/Mu^ = 0, (32) 



or 


(33) 


Analytic solution for the ISCO radius can be found in 
the paper [^. Solutions of (1^ should be substituted 
into (1^ and (1^ for finding E and L. 

In the limit a = 0, we have risco = 6M in the 
Schwarzschild case. In the extreme Kerr limit a = M 
we obtain 


’’ISCO = QM, Lisco = 


22 


M, 


75 ISCO = 


5 


(34) 


for the counter-rotating orbit and 


2 1 

’’ISCO = M, A ISCO = 75ISCO = 

for the co-rotating case 0; 0- 


(35) 


IV. MOTION OF A SPINNING TEST BODY IN 
THE EQUATORIAL PLANE OF A BH 


For a description of the influence of body’s spin on the 
parameters of its circular orbits we use the Mathisson- 
Papapetrou-Dixon (MPD) [IMllil equations of mo¬ 
tion: 


Dt 

Dt 




= pi^v'' -p'^v>^. 


(36) 
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Here D/Dt is a covariant derivative along the parti¬ 
cle trajectory, t is an affine parameter of the orbit 
0 , is the Riemannian tensor, and are 4- 

momentum and 4-velocity of a test body, is its spin- 
tensor. The equations were derived under the assumption 
that characteristic radius of the spinning particle is much 
smaller than the curvature scale of a background space- 
time [13 (see also 0 , 0 ) and the mass of a spinning 
body is much less than that of BH. 

These equations are, however, incomplete, because 
they do not define which point on the test body is used 
for spin and trajectory measurements, so we need some 
extra condition to do that 0. We use the condition of 
Tulczyjew [i3113, that fixes the centre of mass of the 
test body 


= 0. (37) 

Let us consider the motion of a spinning particle in 
the equatorial plane of Kerr BH. In this case the angular 
momentum of a spinning particle is always perpendicular 
to the equatorial plane |24| . Therefore we can describe 
the test particle spin by only one constant s which is the 
specific spin angular momentum of the particle. Value 
|s| indicates the magnitude of the spin and s itself is its 
projection on the z-axis. It is more evident to think of 
the spin in terms of the particle’s spin angular momen¬ 
tum Si = smz which is parallel to the BH spin angular 
momentum S 2 = oMz, when s > 0 , and antiparallel, 
when s < 0. Here z is a unit vector in the direction of 
the z-axis and m is a mass of the particle , [13 . 

Saijo et al 0 have derived the equations of motion of 
a spinning test particle for the equatorial plane of Kerr 
BH. The equations of motion for the variables r, t, tp in 
this case have the form [0 


angular momentum per unit particle rest mass which is 
collinear to the spin of a BH. 

Equation for radial motion can be rewritten as 

(E,A,r)2 = asE^ - 2psE + 7 ^, (40) 


where 


Cts = 


Ps = 


2 2 ^^(r -I- M) 

+ a^ -\ -i- - 


- A(a -I- s)^. 


Ms\ f 2 2 as(r-I-M) 


r + a + 


— A(a -I- s) 


J, 


Ms 


7 s = a -I-- A 


rUl- 


Ms‘ 


+ r 


(41) 


This form is used by Favata [39j- Coefficients as, /3s, 7s 
agree with Rasband 0 results. 


V. THE EFFECTIVE POTENTIAL AND 
CIRCULAR ORBITS FOR SPINNING TEST 
PARTICLES 


The question of what we should consider as the EP in 
the case when a particle has spin does not differ from the 
’’spinless case” in the Kerr space-time. If we want to keep 
analogy with the EP in the Schwarzschild background, 
we should define it as the solution of equation OsE^ — 
2l3sE -I- 7 s = 0 with Os, /3s, 7 s from (gT]). 

EP in this form can be found in the paper 24| : 


Usk{r; J,s) 


+ \//3s - Qs7s 

CXs 


(42) 


(EsAsr)^ = Rs, 

= a f 1 H ——j [J — (a + s)E] H-—— 

EsAs (/3 = ^1 -I—^ [J — (a -f s)E] + "^Ps- 

(38) 


where 


Ss 

As 

Rs 




3Ms^r[—(a -I- s)E + J]^ 

Sf ’ 

p2-A|^ + [-(a + s)i3 + J]2|, 


(39) 


with as, /3s, 7s from (l4l1l. see also Favata 0, Rasband 
171, Suzuki & Maeda0|, 01, Steinhoff & Puetzfeld 
41| . For a review of misprints and comparison of effective 
potentials presented in different works see Appendices of 
the paper [4l|. 

Using this EP for finding the ISCO parameters is pre¬ 
sented, for example, in the paper of Favata 01. It is 
shown that conditions for a circular orbit (the vanishing 
of r and r) are equivalent to 


E = Usk, 


dUsk 
dr 


For the ISCO the third condition is added: 


d^Usk 

dr'^ 


(43) 


(44) 


Ps = 


r'^ + a'^ + 


as(r + M) 
r 


E- 



Here x = dx/dr and the affine parameter r is normalised 
as p^Vi, = —m 01; E is the conserved energy per unit 
particle rest mass, and J = Jz is the conserved total 


A general solution of these three equations for the ISCO 
of a spinning particle is supposed to be found numerically. 
In the small-spin limit, Favata [0 analytically had found 
parameters of ISCO in the Schwarzschild metric. 

Since our purpose is to obtain the analytical expres¬ 
sions for circular orbits in the Kerr metric, it is more 
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convenient to define the ’effective potential’ similar to 
(Ei), by analogy with the work [^. We rewrite equation 
(|40|) as: 


r2 = 


E?A2 


(cXsE"^ — 2PsE + 7s), 


(45) 


and define the effective potential as its right-hand side: 


Vsir;J,E) = 




{asE^ -2PsE+ -/,). (46) 


It can be easily shown that solution of system of equa¬ 
tions determining circular orbits is not changed if the 
effective potential is multiplied by some function which 
is not equal to zero (with its first and second derivatives) 
for physically relevant parameter values. Therefore it is 
more convenient for analytical calculations to define ef¬ 
fective potential as 


V,{r-,J,E) = \iasE^-2PsE + js), (47) 

with as, Ps, Is from (SID- Notice that in the case of 
a spinning particle we use the total angular momentum 
J instead of the orbital angular momentum L0 Further, 
for the sake of convenience, we shall change variables and 
work not with r and J but with u = Ijr and x = J — aE, 
so the function Vs{u;x,E) will be used. 

In order to find the circular orbits we need to solve the 
system of equations 


r ^s = o, 

1^=0 

I du 

The explicit form of these equations is given in the first 
two equations in the system ([50]). The solution of this 
system defines the set of parameters x and E for stable 


and unstable circular orbits as functions of s and u. In 
Figure |2| we present the numerically calculated depen¬ 
dence of r on total angular momentum J and the influ¬ 
ence of spin on it for the Schwarzschild background, see 
also [TtII . 

The point at which upper and lower branches on Figure 
|2| meet defines the ISCO for the given value of spin. In 
order to find the dependence of the radius of the ISCO on 
spin, we need to add the condition d’^Vg/dr'^ = 0, which 
is equivalent to 



FIG. 2. Dependence of the radii of the stable (upper branches) 
and unstable (lower branches) circular orbits on the total an¬ 
gular momentum, for different values of spin s. All parameters 
are given in units of M, see also 0- 


d'^Vg 

dv? 


= 0 , 


(49) 


to the system (l4^ . The resulting system for finding the 
ISCO parameters, thus, consists of three equations (|4^ 
and (1491) and has the explicit form 


' (1 -I- 2asv?' — + 2Ms^u^)E^ + {—2au^x + 2su^x — QMsu^x — 2aMs^u^x)E — 1 -|- 2Mu — + 2Ms^u^ 

- + 2a^Ms^u^ - + 2M^sW - a^M^s'^u^ - u^x^ + 2Mu^x‘^ + 2aMsu^x^ + M'^s'^u^x^ = 0 

(4asM — 2s^u + QMs^v?)E^ {—Aaux + Asux — ISMsu^x — \QaMs^u^x)E + 2M — 2a^u + QMs^v? — I6M^s 
-I- lOa^Ms^'u"' — QM'^s'^u^ + lAM^s'^u^ — — 2ux^ + QMu^x^ + lOaMsu'^x'^ + QM^s^u^x^ = 0, 

(4as — 2s^ -|- V2Ms^u)E^ + (—Aax + Asx — 36Msux — A0aMs^u^x)E — 2(a^ — 6Ms^u + 24M^s^u^ — 20a^Ms 
+ - A2M^s^u^ + 28a^M^s\^ + x^ - QMux^ - 20aMsu^x^ - IbU^s^u^x"^) = 0. 






(50) 


These three equations form a closed system on three pa¬ 
rameters of ISCO E, X and u, which are then dependent 


2 For a spinless particle the conserved quantity is the orbital angu- 
lar momentum whereas in the case of a spinning particle the 
conserved quantity is the total angular momentum J 2 , which in¬ 


cludes spin terms [^ . In this case ’orbital angular momentum’ 
at infinity Lz can also be introduced as Lz = Jz — s, see . 


















only on the Kerr parameter a and particle’s spin s. This 
system can be used for numerical calculation of r, E, J 
of ISCO at given a and s. 

Using EP in the form Usk{r', J) (1^ permits of splitting 
that system of three equations into two parts: equations 
dUsk/dr = 0 and d^Usk/dr'^ = 0 contain only r and J, 
and can be solved independently from equation E = Usk ■ 
So actually the problem is reduced to solving of the sys¬ 
tem of two equations. It is clearly seen in the case of 
the Schwarzschild metric for a spinless particle, see Sec¬ 
tion 3. In the case of using EP in the form Vs{u]x,E) 
6Zl) the equations dVs/du = 0 and d^Vs/du^ = 0 contain 
not only u and cc, but also E, so that in order to find 
these three parameters it is necessary to solve simultane¬ 
ously three equations. Disadvantage of Uskir; J) is that it 
contains radicals, and each differentiation increases com¬ 
plexity of expressions. EP in the form Vs(u]x,E), on 
the contrary, has polynomial form with respect to u, and 
each differentation makes the expressions simpler. All 
these properties are present in the spinless case also, and 
it is convenient to look for the ISCO parameters for the 
Schwarzschild case using Uschwir; L) and the ISCO pa¬ 
rameters for the Kerr case using Vs(u\x,E). 

Another important characteristics of the particle circu¬ 
lar motion is its angular velocity. We calculate the spin 
corrections for the orbital angular frequency of the par¬ 
ticle at the ISCO, as seen from an observer at infinity. 
This angular frequency is defined as 


^ dip/dr 
dt/dr 


(51) 


The values dip/dr and dt/dr are found from the second 
and the third equations in (l38ll . where we should substi¬ 
tute values of r, E and J at a given orbit. To find the 
ISCO frequency ft rsm we need to use the ISCO values 
of r, E and J, see [^ . 


VI. SMALL-SPIN LINEAR CORRECTIONS FOR 

ISCO PARAMETERS AT ARBITRARY KERR 
PARAMETER 

The system (1501) does not admit of an analytical solu¬ 
tion for arbitrary values of a and s. Using s <C M, we 
can expand the equations of the system (EOl) into series 
of powers of s. 

In this section we find analytical solution for ISCO pa¬ 
rameters with small-spin linear corrections, for arbitrary 
value of a. 

To perform analytical calculations and get analytical 
results for small-spin corrections, it is more convenient to 


Further we can say that in both cases total angular momentum 
is conserved with a note that in spinless case it consists of orbital 
angular momentum part only. 


use the following equivalent system of eqnations instead 
of the system (|48ll and ([49]): 


2 du 


= 0 , 


l^u-Vs=0, 

2 du 


Id^Vs 
2 du^ 


= 0 . 


(52) 


Similar approach for spinless case is used in (45l| . 

We shall seek for solution of the system (l5^ in the 
form 


a; = xo -I- sxi, 

E = Eo + sEi, (53) 

It = Mo + SUi- 

Here xq, Eq, uq are the values of x, E, u in the spinless 
case, and Xi, Ei, mi are the linear corrections associated 
with spin. The dimensionalities are [x] = [xq] = [M], 
[xi] = 1, [E] = [Eo] = 1, [El] = [M-i], [m] = [mo] = 
[M ^], [mi] = [M ^], [s] = [Mj. Further we introduce 
these expressions into (1521) and linearise it in spin. Since 
the equations of the system are to be satisfied for arbi¬ 
trary values of spin, we obtain the system of six equations 
for six unknowns xq, Eq, mq, Xi, Ai, mi. 

First, there are three equations for three unknowns xq, 
Eq, Mo, at given a: 

{ M — o^uq — XqMo — 2aEQXoUQ + SMxquI = 0 , 

l-E^- Muo + Mxlul = 0, (54) 
— — x/^ — 2aEoXo + QMxqUo = 0. 


These equations describe the parameters of ISCO for the 
spinless case. Solution of the first two equations in this 
system gives us the values of xq and Eq for an arbitrary 
circular orbit of given inverse radius, see (E51) . Solving 
three equations simultaneously, we can obtain the equa¬ 
tion for Mo, see dSl]). 

The solution of equation (|5^ for mq has a complicated 
form. Therefore it is more convenient to express the six 
unknowns not as explicit functions of a but as the explicit 
functions of a and mq, keeping in mind that uq can be 
found from ea. (|32l) at arbitrary a. Before we proceed we 
need to notice that representation of all unknowns via uq 
given below could be rewritten in a different form using 
(I32|). 

To get Xq and Eq for ISCO we can substitute u = uq 
into the expressions for an arbitrary circular orbit (1281) . 
But it is more convenient to obtain directly the connec¬ 
tions between parameters xq, Eq, uq for ISCO from three 
equations (1511) . It is easy to obtain that: 

Xq = ^0 — ^ ~ ■ (55) 

Whereas Eq has positive value and we can write 


Eq = 



(56) 
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the value of Xg can have both signs. From expression for 
Xq in (1^51) . we see that the second term in numerator is 
always bigger (or equal) than the first term, which means 
that the sign of xg is determined by sign in front of the 
second term in numerator. Therefore xq > 0 stands for 
co-rotating orbits and xo < 0 for counter-rotating orbits. 
We can write: 


Xg 


= T 


1 

V^Ug 


(57) 


The ISCO radius r can be found by using r = 1/u: 

1 Ui , , 

r = rg + sri =- s—^. (64) 

Ug Uq 

We calculate also the ISCO frequency for arbitrary 
value of a, as described at the end of section [V] Using 
dSTl), dMl) and o in m we obtain as 

0 = 0o-bs0i, (65) 


For Eg we can also write: 

^ 6Mug — 3uga^ — 1 {6Mug — — l)-\/3 

° dauQXg 6aug 

(58) 

Three equations for small-spin corrections xi, Ei, ui, 
partially simplified with using the third equation in (l5^ 
and relation XgUg = 1/3, are: 


' (— 6 x 0 ~ 6aEg + ISMxgUg) Xi — 6axg Ei + 

+ GaEg + 6EgXg — 27EgMxgUg + 5aMug = 0 , 

< 2MxgUg Xi — 2Eg El + MuUq — SEgMxgUQ — 0 , 
(—3x0 ~ 3aUo -I- ISMxqUo) Xi — 3axo Ei + 9Mxg ui-|- 
, -I- 3aUg -I- 3EgXg — 27EgMxgUg + lOaMug = 0 . 

(59) 

After significant simplifications with using expressions 
(1551) . (1571) . (1551) we succeed to obtain the corrections in a 
simple form: 


.1 1 
El = —XgMUg = ±—=Mu, 


Xi = — 


Ui = - 


° V3 

(1 — 6Mug + 9a^ttg)xo 
4a 

Ug{l — 6Mug + 5a^Ug) 
2a 


^ 1 — 6Mug + 9a^Ug 


4^, 


'Q,Uq 


(60) 


To avoid singularity in the Schwarzschild case (a = 0, 
Ug = 1/6M), we can also rewrite (1601) as 


El 


— XgMug = E—=Mu 


1 


2 

O’ 


XI = ±—^(3aUo ± 2y/Mug) = —J=(2\/Mug ± 3aug), 
v 3 V 3 


ui = —4ito(auo ± \/ Mug). 

(61) 

Let US summarise the results of this section. To obtain 
ISCO parameters for given a in the form (1531) one should 
find Ug from equation 


1 — 3a^Ug — 6Mug =F 8a Mug = 0, (62) 

and then calculate xq with using (153, Eg with using (1551) , 
xi, El, ui with using dsni) or dm). The value J can be 
found by using J = x -I- aE: 


J — Jg sJi — (xo -l- aEg) s(xi aEi). (63) 


_ \/Muy^ 9'/Mug{\^M ± ay/u^) 

^^0 - - o /9 5 '^1 ~ r-Q ’ 

a^/Mu^/^Tl 2(lTayMaf) 

( 66 ) 

or, in terms of rg as 

^ y/M ^ 9yM(vOT±a) 

“o - 3/2’ -- y— - —(67) 

aVM T rg 2 ^ (rg^^ T a^Mj 

In all the formulae the upper sign corresponds to the 
antiparallel orientation of particle’s angular momentum 
J and BH spin a (counter-rotation, J < 0), the lower - 
to the parallel one (corotation, J > 0); s is the projec¬ 
tion of spin on the z-axis and can be positive (spins of 
particle and BH are parallel) or negative (antiparallel); 
a is positive or equal to zero; the z-axis is chosen to be 
parallel to BH spin a. 


VII. SMALL-SPIN CORRECTIONS FOR ISCO 
PARAMETERS IN SCHWARZSCHILD METRIC 


In order to illustrate the obtained results, let us con¬ 
sider some important limit cases: Schwarzschild metric 
and Kerr metric, when the BH rotation is slow (a <C M) 
and extremely fast (a —)■ M). The latter case is called 
” extreme Kerr BH”, for at the values of a exceeding M 
the BH disappears and we have a naked singularity (T7j | 
without an event horizon. 

In the Schwarzschild case there is no prefential direc¬ 
tion associated with BH spin (a = 0), so z-direction can 
be chosen arbitrarily. Choosing z-axis along J vector, 
and taking ug = IjGM, we obtain: 


<7iSCO 

L^isco 

Cisco 

ISCO 


2'/3M -\— 

2^2 1 sj 

~3 36^3 

6M - 

I Sj 

6V6M 48M2 ■ 


( 68 ) 


Here instead of s, which is the projection on z-axis that 
does not unabiguously correspond to any physical direc¬ 
tion in Schwarzschild case, we use sj which is the pro¬ 
jection of particle’s spin upon the direction of J and is 
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positive when the particle’s spin is parallel to it and neg¬ 
ative when it is antiparallel. Value J is considered as 
positive in this case. 

Small-spin corrections for Schwarzschild metric were 
derived by Favata 


VIII. CORRECTIONS FOR ISCO 
PARAMETERS FOR SLOWLY ROTATING BH 


The case of slowly rotating Kerr BH (a <C M) allows 
us to separate corrections that arise from different types 
of angular momentum coupling. 

Our general formulae, presented in section lyTl allow us 
to obtain correction terms proportional to a, s and as. 
Considering a <C M in equation (1621) . we get ug: 

_ 1 1 [2 a 13 

3M 216M M2’ ' 

Using the first two terms from dMl) in dMl), dSll), dSB), 
we obtain xg, Eg, xi, Ei and ui. Note that when using 
formulae dini) for calculation of corrections Ei and Ui, 
the quadratic terms from should be also taken into 
account. Corrections a, s, as for ISCO parameters can 
be as well obtained directly by solving the system (1501) . 
with keeping terms linear in a and s. 

If we don’t assume any ratio between values a and s, 
a full description of the problem requires also to obtain 
quadratic terms s^ and in ISCO parameters. The 
scheme for obtaining the squared corrections is analogous 
to the linearisation procedure, only we need to express 
the parameters in the form 

X = Xg Sj;i -I- S^X2, 

E = Eg -\- sE\ -\- s^E 2 , (fO) 

U = Ug -\- SUi + S^U2, 


and retain the squared in spin terms in the system (1501) . 

For our purposes we can use more simple way: correc¬ 
tions can be found from solution of the system dSOD 
for Schwarzschild case, and a^-corrections can be found 
from consideration of the spinless case. 

Finally, we obtain: 


J ISCO = t2\/3M- 


± 




27 


± 


r. 2 V 2 
^'isco = —± 


3 
a 

\m) 

1 a 
ISy/lM 


^ ^ 11 a 

3 36^3 M 

1 o 


S± 


I 


± 


5 /a 
162^2 VM 
2 


)■ 


mVlM 

5 


\/2 a s 


432y2M2 ‘ 


I'isco = 6M ± 4 

(s) 


7M 


a±2 
29 


2 a 
9M‘ 


72M 


( 71 ) 


Different situations of mutual orientation of BH and par¬ 
ticle spin are presented on Figure [31 where the case 
a <C M is drawn on the basis of the first three terms 
in formula (1711) for risco- For example, for particle C" 
we should use the upper sign and s > 0, for particle D' 
we should use the lower sign and s < 0. 

This result allows to understand the character o f th e 
spin-orbital coupling in GR. From the terms ±4-\/2/3a 
and ±2^2/35, we see that the the spin-orbital coupling 
depends on mutual orientation of the spin and angular 
momentum vectors: it is attractive and the radius de¬ 
creases, when they are parallel, and repulsive, if they are 
antiparallel. 

The terms with as, a^, s^ in formulae (EID are of second 
infinitesimal order and include information about inter¬ 
actions of particle spin, BH spin and orbital momentum. 
Considering the formula for radius it is easy to show that 
the trinomial containing these three terms is negative for 
all possible values of a and s. Therefore, quadratic cor¬ 
rections always lead to decrease of the ISCO radius. 

Substituting values of JiscOi F'iSCO and risco from 
dZU into dSU, we obtain the ISCO frequency as 

^^isco - + 216MM 48 M 2 ® ^ 

1 as 59 a^ 97 

ISyeM M2 648V6M 3456V6M M2 

(72) 


IX. SMALL-SPIN LINEAR CORRECTIONS FOR 
EXTREME KERR BH 



In the extreme Kerr background (a —>■ M) for a spin¬ 
less particle the ISCO radius is 9M in the antiparallel 
case (counter-rotation), and M in the parallel case (co¬ 
rotation). Numerical solution of the system (l50l) for a 
close to M shows that ISCO radius has very different 
behaviour for co-rotating [J > 0) and counter-rotating 
(J < 0) particles, see Fig. Sj and [5] We will therefore 
consider these cases separately. 

In the antiparallel case we get, using dsiD, dST]), dnoi), 

(ESI): 


Jisco 

FRISCO 

I'lSCO 

^2isco 


22V3,, 82^3 

-M -I- s 

9 243 

5^3 73 s 

~ ^M’ 

16 

9M + -^s, 

1 3s 

~26M 338M2' 


(73) 


For the parallel case we search for ISCO parameters 
at values of a = (1 — d)M, 5 1, because divergences 

may appear at a = 1 at some ways of calculations. Using 
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FIG. 3. Radii of ISCO for different directions of motion and 
spin of test particles. 

Six particles, A, B, C, D, E, F, with different directions of 
spin and orbital motion are presented, for different Kerr pa¬ 
rameter a. For a = 0 we use just letters (e.g.. A), for a M 
we use primes (e.g.. A'), for a = M we use double primes 
(e.g.. A”). Direction of particle’s spin is shown by an arrow 
on a circle: clockwise arrow (C and F) indicates particles 
co-rotating with black hole, anticlockwise arrow {D and E) 
indicates counter-rotating particles, for spinless particles {A 
and B) we use circles without arrow. Direction of orbital mo¬ 
tion around BH is shown by straight arrow outside circles. 
Size of BH is shown for the case of the extreme Kerr space- 
time. The ISCO radii for a <C M are drawn on basis of only 
the first three terms in formula (ED, with s <C a <C M. 

In the case of a Schwarzschild BH (a = 0) the ISCO radius 
for spinless particles equals to 6M, independently of direc¬ 
tion of orbital motion, see A and B. In the case of a Kerr 
BH the ISCO radius depends on direction of orbital motion. 
For slowly rotating BH, a <C M: counter-rotating particle A' 
has orbit with risco > 6M, co-rotating particle B' has or¬ 
bit with risco < 6M. For the extreme Kerr BH, a = M, a 
counter-rotating particle A” has orbit with risco = 9M, a 
co-rotating particle B" has orbit with risco = M. 

Presence of spin increases or decreases ISCO radius in com¬ 
parison with spinless case, see particles C, D, E, F (also 
C', C" and so on). Remind that z-axis is directed along the 
axis of rotation of BH, so a > 0. Total angular momentum 
J = Jz\ J > 0 corresponds to co-rotating of orbit, J < 0 cor¬ 
responds to counter-rotating orbit. Spin s > 0 corresponds 
to spin vector parallel to a (corotation with BH), spin s < 0 
is antiparallel to it. The ISCO radius of spinning particle in 
the Schwarzschild case depends on sj, projection of spin vec¬ 
tor on total angular momentum vector. Therefore particles 
C and D (also E and F), which have different ISCO radii in 
the Kerr metric, have the same ISCO radius in Schwarzschild 
case. For the extreme Kerr BH particles B”, D", F" have 
the same ISCO radii. 


(EH), EH), EQl), E9, we get: 


-^ISCO 


^^isco 


f ISCO 

f^isco 


2 

7!^ 

2 

”7! 

1 

7!^ 

1 

” 7 ! 


Me 


2 X 2^/351/3\ 

7 ! ) 

4 X 2^/3j1/3N^ 

^ TT 

22/3^1/3 

2x2^/3j1/3\ g 

^ 75 ) 


-) 


+ 


(1 + 22/351/3) M - 2 X 22/351/35, 


I 3 X 22/351/3 9 X 22/351/3 

W 8M 16M2 ^ 


(74) 


We see that in the case of a = M (5 = 0) the corrections, 
linear in spin, are absent in formulae for ISCO radius 
and frequency. This was also demonstrated in [l9j |. In 
the next section we discuss this feature in details and find 
the exact in spin expressions for the ISCO parameters. 

The radius of horizon of the Kerr BH, 


Thor = M -b v'M2-a2, (75) 

in the case a=(l — 5)M will have the form ii 

rho, = M(1 -b v^). (76) 

As s <C M and 51/2 < 5i/3 for 5 <C 1, it will always be 
that risco > rhor- 


X. EXACT SOLUTION FOR EXTREME KERR 
BH IN THE CASE OF CO-ROTATION 


In the work on basis of the numerical calculation, it 
was noticed that in the extreme Kerr background for the 
parallel case the magnitude of test-body’s spin does not 
influence the radius of the last stable orbit and it always 
remains equal to M. Solving numerically Eqs EH , we 
plot Fig. |4]where this conclusion is distinctly visible. We 
also succeeded in proving this analytically. 

On basis of results of the numerical solution we have 
substituted u = 1/M into the system (1501) and have ob¬ 
tained the expressions for x and A, which is the analyti¬ 
cal proof of existence of such solution of the system. The 
exact solution is then 


JlSCO 


M2 - s2 

2- , 

My'S-bBs/M 


Eisco 

risco 


M^-s^ 

~ M2y/3-b6s/M’ 
= M. 


(77) 


Using ED in ED, we obtain the ISCO frequency as 




I 

ISCO - 


(78) 
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FIG. 4. The ISCO radius as a function of particle’s spin for 
different values of a close to M: a/M = 0.9 (a), 0.99 (b), 0.999 
(c), 0.9999 (d), 0.99999 (e), in case of co-rotation. In limit 
a —^ M curves tend to the horizontal asymptot risco = 1. 



FIG. 5. The ISGO radii as a function of particle’s spin for 
different values of a close to M: a/M — 0.9 (a), 0.99 (b), 
0.999 (c), in case of counter-rotation. In limit a ^ M the 
ISCO radius at s = 0 tends to the 9M. 


Developing (EZl) in series in s, we get the linear in spin 
terms which are in agreement with (ITU) at <5 = 0. This 
result allows us to obtain the full in spin correction to the 
nearly extreme black hole, with a = {1 — S)M, considered 
in the section El For this purpose after substituting this 
value for a we develop the system (1501) in series in S and 
seek for the solution in the form y = yois) + yi{s)S^^^, 
where y stands for x,E,r. After some algebraic calcu¬ 
lation we obtain for the nearly extreme Kerr metric the 
expressions: 


•Zisco = 2MAisco j 

M2 - s2 


E 


ISCO = 


^isco 


n 


ISCO = 


(m2 - s 2)1/3(2M -b sf/^ ^1/3 

M2y^3-b6s/M y3M5/2(M-b2s)3/2 

, M(M2-s2)1/3(2M + s)2/3 

^ Z{M,sY/^ 

3(M-s)3/3(M-b2s) 


2M 4(2M-bs)i/3(M-bs)2/3Z(M,s)i/3 
where Z{M, s) = + IM^s + 9M^s^ + llMs^ - s'‘. 




(79) 


Developing this result in series in s, one obtains all terms 
in expressions dZil). 

Note that derivation of MPD equations was performed 
under condition s <C M, so always we have Tisco > Oj 
Aisco > 0 in dZZl) and dZSl). Indeed, equations of motion 
of a spinning particle were derived under the assumption 
that its characteristic size is very small compared with 
the characteristic length of the background field (e.g. the 
distance from the central body) [n|, Mpller 

had shown in the context of special relativity that the 
classical spinning body must have a certain minimum 


characteristic size connected with its specific spin, see 
also Wald [S^. For particle at ISCO, in our system of 
units the above conditions reduce to s <C M, see also 

M, M- 


XI. SPIN CORRECTIONS FOR ARBITRARY 
CIRCULAR ORBITS 

In previous parts of the present paper we have dis¬ 
cussed ISCO, and investigated how the presence of spin 
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influences parameters of ISCO for different given a. Let 
us now discuss not innermost stable but arbitrary circu¬ 
lar orbits. We will investigate how the presence of spin 
influences orbit of given radius - namely, how E, J, 
of a circular orbit of given radius are changed if a test 
particle has spin. 

Let us suppose that inverse radius u of circular orbit 
is known and consider the system (l4^ . We shall seek for 

I 


its solution in the form 

X = Xq + SXi, 

E = Eq + sEi. 


(80) 


Here Xq and Eq are the values of x and E in the ’’spin¬ 
less” case, see (l28l) . xi and Ei are the linear corrections 
associated with spin. We obtain for the linear correc¬ 
tions: 


2 EgXo — 9 MEqXou — ^o^MxqU^ -I- 2 aEQ + 8 aME qXqU^ 
2 {aEQ + EqXq — 8 M EqXqu + aM XqU^) ’ 

Mu^xo{—EoXo + ia? EqXqv? — aE^ — AaMxgU^ + SaxQU^) 
2{aE^ + EqXq — SMEqXqU + qMxqU^) 


These expressions are valid for both stable and unstable 
circular orbits with the zero order ISCO parameters given 
by (1^ . Value of total angular momentum J = Jq + sJi 


can be obtained using J = x + aE. 

The expressions (ISTl) can be simplified further. Using 
(E51) we obtain: 


Mu^/'^{ay/u ± '/M){1 -\- ± Aq'/MtE) 

2(1 - 3Mu T 2a\/Mu3)3/2 


(82) 


J _ SMa'^u^ ± -k T -k 2Mv?c? ± 21M^/^u^/‘^a + T 9ay/M' 

2(1 - 3Mm =f 2aVMu3)3/2 


m3 _ 13Mm -k 2 


(83) 


Written as functions of r the resulting formulae are 
^3/2 _ =F 

° ^ r3/4(r3/2 _ 3Mri/2 ^ 2aMi/2)i/2 ’ 

± -k a^) 

= r3/4(r3/2 _ 3Mri/2 ^ ’ 


J 


Ui 




M(a ± \/Mr){r^ -k 3a^ ± Aa^fMr) 

2rii/4(r3/2 _ 3M=F 2ayM)3/2 ’ 

2r^ — 13Mr‘* =k + 18M^r^ ± 2\aM^/'^r^/'^ -k 2a^Mr^ =k + Aa^M'^r ± -k 


2rii/4(r3/2 - 3M=F 2a\/M)3/2 


(85) 


Substituting (|M|l and (IM)l into ([5T|l . we obtain angular 
frequency U = Uo + s^i at the circular orbit with a given 


inverse radius u as 


VMu^ 
aVWv^ =F 1 

8u^'/M{'/M ± Uy/u) 
2{aVMu^ =F 1)^ 


(86) 
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or, as function of radial coordinate r, as 


fin = 


Vm 


= — 


U'/M =F r3/2 ’ 

± a) 


2i/r(a-\/M =F j-3/2^2 


(87) 


The linear corrections for the energy, angular momen¬ 
tum and frequency at a g iven radius for the circular or¬ 
bits are given also in jl9j |. in another presentation, 
see also Appendices A and B in [4l| where comparison of 
spin orientation in different papers is discussed. 


or a test-body are parallel and repulsive when they are 
antiparallel. 

(iv) Linear in spin corrections to the ISCO parameters 
for extreme Kerr BH both for counter- and co-rotating 
cases are obtained, see (ES]) and dZil). 

(v) The exact parameters of ISCO for a spinning test- 
body in the extreme Kerr background in corotation case 
are found, see (EZl) and (17^ . It is proved analytically 
that radius and angular frequency of such orbits are in¬ 
dependent of the particle’s spin while the values of energy 
and total angular momentum depend on it. 

(vi) Figures [T] and [3] illustrate the number of possible 
configuration for rotating BH and a spinning test particle 
with parameters calculated in the present work. 


XII. CONCLUSIONS 

(i) The linear in spin corrections for the ISCO param¬ 
eters: radius, total angular momentum, energy, orbital 
angular frequency for arbitrary Kerr parameter a, are 
found analytically, see Section IVll 

(ii) Small-spin linear corrections for Schwarzschild 
metric (firstly obtained by Favata [s^) were derived, see 
formulae (l68l) . 

(hi) Expressions for the ISCO parameters of a slowly- 
rotating Kerr BH up to second order in spin and Kerr 
parameter are obtained, see (El). The case of a slowly 
rotating Kerr BH gives the idea of spin-orbital coupling’s 
influence on the parameters of ISCO. It is found that 
spin-orbital coupling has an attractive character when 
the total angular momentum J and spin of either BH 
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